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I. Introduction 

To reduce je,t-englne-produced NOx and CO sisiultaneously, 
one of the advanced combustor concepts Invloves lean premlxed/ 
preyaporired combustion system. In the systet;>, since liquid 
fuel droplets evaporation in present-day combustion chambers 
is absent, locally fuel-rich and high temperature regions are 
removed. As a result, substantial reduction of NOx is possible 
without. unacceptable CO emissions. Furthermore, because of 
lower flame temperature and more uniform combustion , smoke and 
radiation may be reduced, turbine inlet temperature may be 
ii!^roved and turbine erosion from carbon particles may be 
avoided . ^ 

Acompanying lean premixed/prevaporized combustion system, 
however, are problems such as much narrower flame stabilization 
limit, which result in blowoff, ignition and altitude relicrht 
difficulties, and flash-back and auto-ignition, which result in 
the presence of flame in the premixino/prevaporizinq section 
with potential subsequent hardware damage. 

To understand the flame flash-back phenomena more extensively, 
an analytical model has been formed and a numerical program has 
been written and tested to solve the set of differential 
equations describing the model. Initial results show that under 
a given set of conditions flame propagates in the boundary layer 
on a flat plate when the free stream is at or below 1.8 m/s. 

More effort is needed in working out with a turbulant reaction 
model and conducting various parametric studies. 
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The formulation of the problem, a general discussion of 
the finite difference solution of the governing equations and 
the results so far obtained and their discussion are contained 
in the following sections. 


II. Formulation of the Problem 


As shown in Figure 1, it is assvimed that a uniform flow 
proceeds from left to right passing on top of a flat plate. 
Downstream a flame is established by a hot segment of the plate. 
Cartesian coordinates are used for the flow field. 

A. Equations of Continuity and Momentxim 
Continuity 

Momentum 


<>t 

<?t 






In the above, u and are x- and y-direction component 
velocities respectively , is the effective viscosity which is 
discussed in II. E. below, J the density, t the time, and p 
the pressure. 

To obtain the variation of ^ , the equation of the state 


is used. 




tfhere i»F (fuel) or 0 (oxygen) ; Yi Is the species mass function; 


^ the density, f Is the effective transport coefficient; D 
is the diffusion coefficient, and W1 is the tintt-mean volumetric 
rate of fuel consun^tlon depending on the combustion model 
described in section II. F. below. 


C. Enthalpy Equation 

Where is the heat of combustion; Cp specific heat of the 
gas mixture; T the temperature and 7\ the conductivity; and 


h * + Hfi( Yp 

D. Boundary Conditions 

Boundary conditions are applied at all form boundaries as 
described below. 

Far upstream, u is specified, !>■ is zero, the fuel fraction is 
0.059, and the enthalpy is calculated from the specified heat of 
combustion and ambient temperature. The incoming air is 
stoichiometrically mixed with methane. 

At the upper boundary, u, p and h have Nexunann conditions, 

V is zero and Yp has Dirichlet condition there. 

At the lower boundary u has Neumann and Dirichlet conditions 
before and after the leading edge of the flat plate respectively, 
b is zero and p, h and Yp have Neumann conditions except that 
along the hot segment of the plate h has Dirichlet condition 
there . 







At the downstream end, V Is zero, p, h and Yp have Neumann 
conditions, and u is modified such that the mass flow rate is 
consistent with its incoming rate far upstream, 

E. Transport Coefficients 

It is assumed that the transport of momentum, mass and 
heat by molecular and turbulent mixing obeys laws of the 
"coeff icient-times-gradient" type associated with the names 
of Newton, Picks and Fourier. The subscripts "eff" implies 
that the turbulent fluid is being treated as a laminar one 
with a set of transport coefficients having augmented values. 
Strictly speaking, the "effective" transport coefficients may 
be expected to posses directional properties: for turbulence 

is rarely isotropic in practically interesting situations. 

However physical knowledge of these matters is still at an 
early stage. We shall assume isotropicity as frequently treated, 
a) The quantity effective viscosity is defined as 

where the eddy viscosity for momentum, assumes the 
well-known Prandtl mixing-length model, i.e. 



where is the mixina length and, in the viscous region, it 
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is expressed, according to Van Driest , by 
j = k ^ [ I - J 


Xn th« above «a^r«Mionf ^ ^ 

71 ss where the subscript o aeana 

a 

that the quantity Is evaluated at th«t wall, K ■ 0.41, and 
A*** ■ 25 (Reference No. 8) . 

b) To express the quantity, one defines a turbulent 

heat conduction in parallel with Fourier’s law of heat 
conduction 

Then one may form an effective heat conduction as 
By defining a turbulent Prandtl number as 

&t= 


it can be shown that 



JL 

B- 



The above equation governs the heat transport coefficients 
in the energy balance equation. = 0.7 and P^^ = 1.0 

Q Q 

are frequently assumed '. 

c) To express the quantity , it is assumed that, for 
both laminar and turbulent flows, Lewis number is eaual 
to 1 as frequently the practice^^, i.e. 

It = ' 

Then, 

= Ip, -Ip, 
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P. Reaction ^todel8 

A one step forward overall gas phase reaction of second 
order is taken to represent the chemistry: 

C^H0 5O2 5aH2 >-^41120 + 3C02 5aH2 where a is the 
molar ratio. 

(a) Laminar reaction rate 

To complete the species equation, 

where F is the pre-exponential constant: p, pressure; 

Yp& Yq, the species concentrations of fuel and oxidizer; 

E, activation energy; and R, the gas constant. 

In the initial results discussed in Section III only 
laminar reaction rate is used. 

(b) Turbulent reaction rate 

The eddy-break-up model^^ for the turbulent flame is 
used as follows. 

First, the reactedness of the gas, "C , is defined as 

7 - " hn _ T- Tu 

^ % -Th 

Then the time-mean volumetric rate of fuel consumption 

may be expressed as 

/ ' « iiil. ■' 

where 


*7i 


-«>- * air fuel mass ratio 
s flame t^nperature 
p B pressure 
c » 0.35 

In this model, the local reaction rate is taken to 
depend not only on a blmolecular Arrenious expression but 
also on the rate of break-up of the eddies by the action 
of turbulence. 


G. Viscosity 

A simple power law^ is used. 

a 

Where yUc - h U53 ^ lO'^ 

Tc 2.73- S \ 


III. Finite Difference Solution of the Governing Equations 


A. Coordinates and Grid Distribution 

A rectangular field of 4.374cm x 1.413cm is imposed on the 
main region of the flow field shown in Fia. 1 . In both x- and 
y-directions, uniform fine grid lines prevail over the interested 
region centered on the established flame front. Far away from 
the flame front, i.e. outside of the interest4e region, the 
grid lines follow a geometric series withthe exception of the 
first and second grid lines near the boundary. 

The magnitude of grid space (space between grid lines) is 
crucial to the proper propagation of the flame. It should be 
a fraction of the flame thickness which can be estimated by 
the thermal diffusivity divided by the reference flow velocity. 


tSt 


In the program, ioc la choaen to be 0.1 mm in the Intereated 
region reaultlng in a grid field of 93 x 31 pointa. 

B. Control Volume 

Figure 2 ahows a grid point and ita control volume. The 
arrangement ia baaed on the work of Patankar and Spalding. 

The control volume, indicated by the dashed lines, is the 
region traversed by a rectangle of sides ^ and The area 

of the control surface normal to the x-direction is simply 
since a unit thickness is ass\uned everywhere in the two dimensional 
region. The surface normal to the X"direction and located at e 
is midway between p and E. As a result p is not always the mid- 
point of e and w. 

With the above arrangroent a staggered grid which is similar 
to the one used by Harlow et. al.^^ is then used. It has the 
convenient feature that the velocity components, a and l> , are 
stored at just the points at which they are needed for the cal- 
culation of the convective contribution to the balance of the 
quantity of question. 

A finite difference equation for a node P is calculated by 
integrating the differential equation over the control volume 
indicated in the above. For evaluating the convection and 
diffusion fluxes through a control volume surface the arithmetic 
mean of the property values on either side of that face is taken. 
For other purposes a step-wise variation with discontinuties at 
the control volume surfaces is assumed. 
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C. Modifications to Finite Difference Equations 

To ensure that the difference equations behave realistically 
when the particular method of computation leads to physically 
unreasonable results, the high-lateral-flux modification^^ is 
used. Under-relaxation is also used to ensure proper convergence 
of the solution. 


D. Solution of the Finite-Difference Equations by TDMA 

The finite-difference equations are solved by the successive 
use of the tri-diagonal matrix algorithm in the x and y directions. 
The TDMA used is similar to the one developed by Roache.^^ The 
interior-point equation is written as 


The solution is of the following form 
nJ; = -r f3; 


^ I y H 


where 


A.' 




and 


C: 6^. g. 

p , — B: A',-i 


For the left-hand boundary condition, one use; 
(a) Dirichlet condition 



(b) Neumann condition 



For the right-hand boundary condition , 


(a) Dirichlat condition 

(b) Nainaann condition 
«M ■ 

Results and Discussion 

The values of various parameters in the numerical confutation 
are listed in Table 1. The programs were run on IBM 3033 at 
Princeton University Computer Center. 

It is assiuned that the steady stream on the flat plate in 
Figure 1 is a stoichiometric methane - air mixture under 
atmospheric pressure and tenferature. Using Prandtl mixing length 
for the eddy diffusivity for momentum end using the bi-molecular 
Arrhenius expression for the chemical reaction the program was 
run for various upstream velocities to determine approximately 
the maximum flow velocity beyond which a flame, which is initiated 
by the hot segment (at 2000 *K) on the plate, would not propapate 
upstream. It was found that at a critical upstream speed of 
1.8 m/s, the flame propagated upstream and that at 1.9 m/s the 
flame remains stationary. This compares with the experimental 
result of Gross^ who reported that flash-backs occu-red around a 
flow speed of 3 m/s using propane-air gas mixture over an airfoil 
plate. The critical upstream speed in the present program is 
expected to increase if eddy-break-up reaction model is used. 

Figure 3 shows the flame envelope and typical velocity, 
temperature and concentration profiles at x 2.2 cm for an 
upstre 2 un speed of 1.9 m/s. At the entrance (far upstream) it is 


tilt 


at one atmospheric pressure of 101.325 while at the exit 

2 

(Far down-stream) 101.189 KN/m . Near the upper boundary the 
fuel CO icentration approaches ts the up-stream value and the 
velocity exceeds the up-stream value by about 10%. The temp- 
erature profile shows a maximtim in the middle but returns to tae 
fixed wall temperature of 2000^K on the plate on the one side 
and approaches the up-stream temperature on the other. 

To reach a steady state from a given set of initial profiles 
such as the one in the above 40 time steps ( t^O.l sec) and 1 
minute and 35 seconds of computer time 2 ire needed to run the 
program. 

The eddy-break-up reaction model has yet to be incorporated 
into the program and a series of paramatric studies need to be 
conducted. 


RBFElOmCSS 


1. Vmxkanp, F.J.» Vsrdoulv/ A.J. and Tonllnson, J.G., **Xra>act of 

lla 9 ulatioi.'s on Putura Gaa Turbina Engine Condsuators" p. 340, 

J. Aircraft, June 1974. 

2. Narak, C.J., Papathakoa, L.C. and Varbulacz, P.I9., "Preliminary 

Studiaa of Autoignition and Flash-back in a Pramixing/Prevaporizing 
Flana Tube Using Jet-A Fuel at Lean Equivalence Ratios" NASA 
TH X-3526, May 1977. 

3. Blazowski, W.S. and Bresowar, G.E., AFAPL-TR-74-32, May 1974. 

4. Anderson, D.M., NASA TM X-3301, 1975. 

5. Gross, R.A., "Boundary Layer Flame Stabilization" p. 288, Jet 

Propulsion, 1955. 

6. Plee, S.L. and Mellor, A.M., "Preliminary Investigation of Flash- 

back in Prevaporlzlng/Premixing Combustors" HSS/CI 77-22. 

7. Van Driest, E.R., "On Turbulent Flow Near a Wall” p. 1007, J. 

Aeronautical Sciences, November 1956. 

8. Crawford, M.E. and Kays, W.M., NASA Report CR-2742, Dece.'nber 1976. 

9. Schlichting, H. "Boundary 'layer Theory" p. 260, 6th ed. McGraw 

Hill, 1968. 

10. Bird, Steward and Liohtfact, "Transport Phenomena” p. 629, Wiley, 

1960 ed. 

11. Spalding, D.B., "Mixing and Chemical Reaction in Steady Confined 

Turbulent Flames" p. 649, Thirteenth Sym. on Combustion, 1970. 

12. Pantakar, S.V. and Spalding, D.B., "A Calculation Procedure for 

Heat, Mass and Momentum Transfer in Three-dimensional Parabolic 
Flows" Int. J. Heat Mass Transfer, 15, p. 1787-1806, 1972. 

13. Harlow, F.H. and Welch, J.E., "Numerical Calculation of Time- 

dependent Viscous Incompressible Flow of Fluid With Free Surface" 
Physics Fluids, 8, p. 2182, 1965. 

14. Patankar, S.V. and Spalding, D.B., "Heat and Mass Transfer in 

Boundary Layers" 2nd ed. Intertext Books, London, 1970. 

15. Roache, P.J., "Con^utional Fluid Dynamics” Hermosa Publishers, 1976. 


13 


Table 1 

Niunerlcal Values Used In the Program 


Parameter 

A-t- (Defined in Section II. E. (a)) 

Cp 

E Activation energy 
K (Defined in Section II. E. (a)) 

Le Lewis no. 

Pr Prandtl no . 

Prt Prandtl no. for turbulent flow 
Q Heat of combustion 
R Gas constant for air 
Tg Temperature of burned gas 
Ty Temperature of unburned gas 
Yf yFuel concentration in the upstream 
Z Pre-exponential constant 


Numerical Value 
25 

265 cal/Kg-OK 
40000 cal/mole 
0.41 
1 

0.7 

1.0 

l.ioezxio”^ cal/Kg 

286.8 N-M/Kg-OR 

2000°K 

300°K 

0.0594 

0.5 
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